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ABSTRACT 

Recent results on physics in the charm energy region are reviewed. Theo- 
retical puzzles related to the exclusive hadronic decays of J/ ip and ip' are 
analyzed. New results and issues on possible glueball candidates such as 
. ^(2230), ?7(1440) and /o(1500) observed in J /ip radiative decays and other 

^ \ experiments are emphasized. Problems in charmonium production with 

Q^ \ large transverse momentum at the Tevatron and fragmentation mecha- 

^ ■ nisms of quarks and gluons, and in particular the ip' surplus and color- 

Q\ ■ octet fragmentation are discussed. Some theoretical and experimental 

^ . results in charmonium and open charm hadrons are also reported. 
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I. Introduction 



Physics in the charm energy region is in the boundary domain between perturbative 
and nonperturbative QCD. The study of charmonium physics has recently received 
5^ I renewed interest. The observed hadronic decays of charmonium may give new chal- 
lenges to the present theoretical understanding of the decay mechanisms. More glue- 
ball candidates are observed in charmonium radiative decays, and are arousing new 
studies of glueball physics. The observed prompt production of charmonium at the 
Tevatron and the serious disagreement between expected and measured production 
cross sections have led to new theoretical speculations about charmonium spectrum 
and novel production mechanisms. There are also many new results in open charm 
physics, including new measurements of charmed meson and baryon decays. In this 
report some of the new results and the status of theoretical studies of physics in the 
charm energy region will be reviewed. 

II. Problems in Charmonium Hadronic Decays 

Charmonium hadronic decays may provide useful information on understanding the 
nature of quark-gluon interactions and decay mechanisms. They are essentially re- 
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lated to both perturbative and nonperturbative QCD. The mechanism for exclusive 
hadronic decays is still poorly understood. One of the striking observations is the so- 
called "pvr" puzzle, i.e., in the decays of J /ip and ip' into pvr and K*K the branching 
ratios of ip' are greatly suppressed relative to that of J/ip New data from BES not 
only confirmed this observation but also found some new suppressions in the 00/2 and 
pa2 channels [H . This gives a new challenge to the theory of charmonium hadronic 
decays. 

Because for any exclusive hadronic channel h the decay proceeds via the wave 
function at the origin of the cc bound state, one may expect 

~ B{J/i, ^h)^ B{J/^ 3g) ^ B{J/i, -> e+e-) ^ ^ ' 

Most channels like 3(7r+7r-)7r°, 2(7r+7r-)7r^, {'K^'K^)pp, Tc^pp, 2(7r+7r-), and the newly 
measured pp, AA, E'^S'^, by BES seem to approximately respect this relation. 

But Qh for pvr and K*K were found (and confirmed by BES) to be smaller by more 
than an order of magnitude than the normal value 0.14. The new BES data give 

ii 

Qp^ < 0.0028, Qk'±k^ < 0.0048. (2) 

This puzzle has led to some theoretical speculations. 

The Hadron Helicity Conservation theorem in QCD |Q, suggested by Brodsky 
and Lepage, indicates that bacause vector-gluon coupling conserves quark helicity for 
massless quarks, and each hadron's helicity is the sum of the helicity of its valence 
quarks, in hard process the total hadronic helicity is conserved (up to corrections of 
order m/Q or higher) 

E = E ^H. (3) 

initial final 

According to this theorem the decays of J/ ip and tp' — > VP (vector and pseudoscalar 
such as pn and K*K) are forbidden. This seems to be true for %p' but not for J/ip. 
The anomalous aboundance of VP states in J / ip decay then needs an explaination. 

In this connection one possible solution to the pvr puzzle is the J/ip — O mixing 
models These models, though slightly different from each other, have the 

same essence that the enhancement of J/ip ^ pvr, K*K is due to O ^ pvr, K*K, 
where O could be a Pomeron daughter 0] or a vector glueball , which could lie 
in the region close to J/ip mass and then mixed with J/ip but not ip'. 

It has been suggested to search for this vector glueball in processes J/ip, ip' — * 
{rj, T]', vrvr) + O, followed by (9 — > pvr, K*K. Obviously, the J/ip — O mixing model 
depends heavily on the existence of a vector glueball near the J/ip. It is therefore 
crucial to search for it in the vicinity of J/ip. But so far there seem no signs for it. 
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Another proposed solution to this puzzle is the so-called generalized hindered Ml 
transition model It is argued that because J jip ^ rjc'^ is an allowed Ml transition 
while — > 77c7 is hindered (in the nonrelativistic limit), using the vector-dominance 
model to relate ip' — *• 777c to ip' — > iprjc one could find the coupling G^/^r^c is much 
smaller than G^^^ic^ ^ind then by analogy, the coupling Gi^ip.,^ would be much smaller 
than Gi^pT^. Assuming ip' —>■ pvr to proceed via ip' — uj' mixing, while ip p-K via 
'p — uj mixing, one would find that 'p' pvr is much more severely suppressed than 
-p — s> pvr. 

There is another model in which a hadronic form factor is introduced to ex- 
ponentially decrease the two meson decays of -p' relative to J/'p. But this model 
predicts a large suppression for many two meson modes, which may not be compat- 
ible with the present data. There is also a proposal to explain this puzzle based on 



the mechanism of sequential quark pair creation |]T0 . 

Now the new BES data give a new challenge to these speculations. It is found 
that in addition to pvr and K*K the suppression also exists in the VT (vector-tensor) 
channels of 'p' decays such as 'p' i^/2(1270) 00 

Q^f, < 0.022, (4) 

and the preliminary data on pa2, K*K2, 0/2 channels seem to also show suppressions 
for 'p', whereas in the 6^7r^ channel there is no suppression is observed for ■?/''. 

The VT decays do not violate helicity conservation, therefore the suppression is 
hard to understand. Moreover, in the J/'p — O mixing model the O VT decay is 
not expected to be a dominant mode, and therefore J /ip ^ VT may not be enhanced. 
Moreover, using the vector dominance model one might relate ip' —>■ ujf2 to ip' —>■ 7/2, 
but the observed 'p>' — ^ 7/2 is not suppressed ||11|| , and this is also confirmed by BES. 



In the generalized hindered Ml transition model, the coupling G^^i^f^ for uj' ujf2 
should not be suppressed because by analogy the coupling G^'^^^^ is not small due 
to the fact that the El transition -p)' 7Xc2 is not hindered. Therefore via -p)' — u' 
mixing the ip' ^ u' ^ 00/2 decay is expected to be not suppressed. It seems that 
within the scope of proposed models and speculations the puzzles related to the VP 
and VT suppressions have not been satisfactorily solved yet. 

In order to understand the nature of these puzzles, systematical studies on J/ tp 
and -pj' exclusive hadronic decays are needed. Many different decay channels such 
as VP (pTT, K*K, ur], ur]', (prj, (pr]', and isospin violated cjtt", pr], pr]', 07r°,---), 
VT iujf2, pa2, 0/2, Pfl ■ ■ ■) AP(6i7r, ■ ■ ■), TP (a27r, ■ ■ ■), VS(cu/o, 0/o, ■ ■ ■)_, VA (0/i, 
ujfi, ■ ■ ■) and three-body or many-body mesonic decays {utt'^it', p>KK,---) and 
baryonic decays {pp, nn, AA, SS, SS, ■ ■ ■) are worth studying and may be helpful 
to reveal the essence of the puzzle and the nature of decay machnisms. In addition, 
to test the hadron helicity conservation theorem, measurements of the decay angular 
momentum distribution are also important. E.g., it predicts a sin^ 6' distribution for 
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Since the r^c, systems are the counterparts of J/'?/', '?/'', it has been suggested to 



study exclusive hadronic decays of r^^ and 77^ |]T3[|T^. It is argued that for any normal 



hadronic channel /i, based on the same argument as for J/?/' and ^Z'', the following 



relation should hold |T^ 



This relation differs from the "0.14" rule for J/ip, because rj'^ 2g is the over- 
whelmingly dominant decay mode, whereas for ip' the ip' — ^ J/ipTTir and ip' 
IXcj {J = 0,1,2) transitions are dominant. As the "0.14" rule for J/ip and ip' , 
this relation for rjc and t]'^ may serve as a criterion to determine whether there exsit 
anomalous suppressions in the t^c, v'c systems. As pointed out in [|13| that since the 
observed rjc VV{pp, K*K*, (pcj)) and pp decays, which are forbidden by helicity 
conservation, seem to be not suppressed, there might be a trigluonium compo- 
nent mixed in the rjc. It then predicts a severe suppression for these decays of r/^, 
which is not close to and therefore not mixed with the 0^+ trigluonium. The t]c 
and hadronic decays are being searched for at BES/BEPC, and will be studied at 
the r-charm factory in the future. In this connection, it might be interesting to see 
whether E760-E835 experiment can find t]'^ in pp rj'^ 27. If rj'^ — > pp is severely 
suppressed by helicity conservation, as the counterpart of tp' —>■ pvr, then it would 
be hopeless to see rj'^ in pp annihilation. Therefore the E760-E835 experiment will 
further test helicity conservation and shed light on the extented "pvr" puzzle. 

On the other hand, the theoretial understanding for these puzzles and, in gen- 
eral, for the nature of exclusive hadronic decay mechanisms is still very limited. It 
concerns how the cc pair convert into gluons and light quarks and, more importantly, 
how the gluons and quarks hadronize into light hadrons. The hadronization must 
involve long distance effects and is governed by nonperturbative dynamics. These 
problems certainly deserve a thorough investigation in terms of both perturbative 
and nonperturbative QCD. 

III. Search for Glueballs in Charmonium Decays 

Existence of the non-Abelian gluon field is the key hypothesis of QCD, and observa- 
tion of glueballs will be the most direct confirmation of the existence of gluon field. 
Charmonium radiative decays into light hadrons proceed via cc — > 7 + (7 + (7 and 
are then the gluon-rich channels. Therefore, charmonium especially J/ip radiative 
decays are very important processes in the search for glueballs. Recent experimental 
studies indicate that there are at least three possible candidates of glueballs which 
are related to J/ip radiative decays. 
• ^(2230) = (?)++. 



The new data from BES HI 15 confirmed the Mark III result 16 and found four 



decay modes of ^ ^ tt+tt , K^K , KsKs, pp in J/ip 7,^ with a narrow width of 
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^ 20 MeV. The branching ratios are found to he B{J/iIj ^ 'y^) x X) 
(5.6, 3.3, 2.7, 1.5) x 10"^ respectively for X = tt+tt-, K+R-, KsKs,pp. 



Combining these data with the PS 185 experiment on pp ,^(2230) — ^ KK [|T^: 
^ pp) X — »• KK) < 1.5 X 10~^ (for J=2) reveals some distinct features of 
the ^(2230): the very narrow partial decay widths to vrvr and KK (less than 1 MeV 
with banching ratios less than 5%); the large production rate in J/ip radiative decays 
{B{J/ip — 7^ > 2 X 10^'^); the flavor-symmetric couplings to tttc and KK. These 
features make ^(2230) unlikely to be a qq meson but likely to be a J^'" = {even)'^~^ 
glueball [IHIli. 

The ^(2230) once was interpreted as an ss meson[^. But a recent quark model 
calculation ||2l| for decays of 1^F2 and l^F^ ss mesons shows that the widths of 1^F2 
and 1^F4 ss mesons are larger than 400 MeV and 130 MeV respectively. The partial 
width of 1^F4 to KK is predicted to be (14 — 118)MeV, also much larger than that 
of ^(2230). Moreover, the lattice study of SU{3) glueballs by the UKQCD group 
suggests the mass of 2+"*" glueball be 2270 ± lOOMeVj^, consistent with the mass of 
^(2230). But the spin of ^(2230) has not been determined yet. (J^^ = 4++ will not 
favor a glueball because it would require a non-S wave orbital angular momentum 
between the constituent gluons, and then lead to higher mass and lower production 
rate in J/ip radiative decays than ,^(2230)). Moreover, in order to see through the 
nature of .^(2230) (e.g., by further examining the flavor-symmetric decays and the 
difference between glueball and qqg hybrid), more data are needed for other decay 
modes, such as rjrj, rjr]', r]'r]' and TnTirir, irirKK, pp, K*K*, uj(j), (j)(f), etc. 

• ?7(1440) J^^ = 0-+. 

For years this state has been regarded as a good candidate for the O""*" glueball. 
However, since both Mark 111 |^ and DM2 find three structures (two ^ and one 
1++ ) in the energy region 1400- 1500MeV in J/ip -/KK-n-, the status of i/?7(1440) 
as a ^ glueball is somewhat shaky. But the new (preliminary) generalized moment 
analysis of BES ||2^, which avoids the complicated coupling effects from different 
intermediate states {K*K and aovr), indicates that the 7^(1440), being one of the 
three structures, may have a larger production rate in J/ip radiative decays with 
B{J/iIj 77/(1440)) ■ 5(7/(1440) KKn) ^ 2 x IQ-^. This may reinforce the 
7/(1440) being a 0""*" glueball candidate. 

While more data and analyses are needed to clarify the discrepancies between 
Mark III, DM2, and BES, some theoretiacl arguments support t/7/(1440) being a ^ 
glueball. The helicity conservation argument favors ^ glueball decaying predomi- 
nantly to KKtt p6|. Working to lowest order in 1/Nc and using chiral lagrangians 



also get the same conclusion ^Tj. However, the lattice QCD calculation by UKQCD 
predicts the mass of 0~"^ glueball to be ~ 2300MeV |2^, much higher than 1440 
MeV. 

• /o(1500) J^^ = 0++. 

The Crystal Barrel (CBAR) Collaboration El at LEAR has found /o(1500) in 
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pp — > 7r°/o(1500) followed by /o(1500) 7r°7r°, 7777, 7777'. This state might be the 
same particle as that found by WA91 Collaboration in central production pp 
Pf{2'K^2TT~)ps and that found by GAMS in tt~p rjrj'n, rjrin, An'^n, namely, 
the (j(1590) |3^. So far no signals have been seen in J/ip radiative decays in channels 
like TTTT, ?7?7, 7]7]' for /o(1500). However, it is reported recently that re-analysis of Mark 
III data on J/tp — > 7(47r) reveals a resonance with J^^ = 0++ at 1505 MeV, which has 
a strong aa decay mode ||3l|. If this result is confirmed, /o(1500) may have been seen 
in three gluon-rich processes, i.e., the pp annihilation, the central production with 
double Pomeron exchange, and the J/ip radiative decays, and is therefore a good 
candidate for 0"*"^ glueball. It will be interesting to see whether /o(1500) aa —>■ An 
is the main decay mode in the CBAR experiment. The mass of /o(1500) is consistent 



with the UKQCD lattice calculation |22 



The theoretical understanding of glueballs is still rather limited. There are differ- 
ent arguments regarding whether glueball decays are flavor-symmetric. 
• Helicity Conservention. 



It was argued by Chanowitz p6| that although glueballs are SU{3) flavor singlets 
it is inadequate to use this as a criterion for identifying them because large SU{3) 
breaking may affect their decays. In lowest order perturbation theory the decay 
amplitude is expected to be proportional to the quark mass 

M{gg qq)j=o oc m,, (6) 

so that decays to ss are much stronger than uu + dd for 0+"'" and O""*" glueballs. This is 
a consequence of "helicity conservation"- the same reason that r(7r ^ fiu) ^ r(7r ^ 
eu), and this might explain why L/ri{1440) — KKn is dominant. 
• Discoloring of gluons by gluons. 



It was argued by Gershtein et al. that due to QCD axial anomaly the matrix 
element as< 0\GG\ri' > gets a large value. Therefore, if the glueball decay proceeds 
via production of a pair of gluons from the vacuum and recombination of the gluons 
in the initial state with the produced gluons, then decays into t]' will be favored. This 
may explain why the O^"*" (^(1590) has a larger decay rate into rjr]' than rji], tttt, KK. 

• Glueball-gg mixing. 

It was argued by Amsler and Close ||3^ that for a pure glueball Go flavor democ- 
racy (equal gluon couplings to uu, dd and ss) will lead to the relative decay branching 
ratios tttt : KK : r]r] : rji]' = 3:4:1:0. Then by mixing with nearby qq isoscalars 
the mixed glueball state becomes 

\G >= \Go > +^{\uu > +\dd > +uj\ss >), (7) 

and the observed decay branching ratios tttt : KK : rji] : 7777' = 1 :< 0.1 : 0.27 : 0.19 for 
/o(1500) may be explained in a color flux-tube model with certain values for mixing 
angles ^ and u with the nearby /o(1370) (an uu + dd state) and an ss state in the 
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1600 MeV region. The problem for /o(1500) and /o(1710) has also been discussed in 
Ref.fBl. 



• Resemblance to charmonium decays. 

It was argued |T8| that pure glueball decays may bear resemblance to charmonium 
decays, e.g., to the Xco{0~^~^) and Xc2{'2~^~^) decays. Both Xco and Xc2 decays may 
proceed via two steps: first the cc pair annihilate into two gluons, and then the two 
gluons hadronize into light mesons and baryons. The gluon hadronization appears 
to be flavor-symmetric. This is supported by the Xco and Xc2 decays, e.g., Xco is 
found to have the same decay rate to tt+tt" as to K~^K~, and the same decay rate 
to 7r"'"7r~7r"'"7r~ as to 'n'~^7i~K~^K~ , and this is also true for Xc2 decays. For a glueball, 
say, a 2++ glueball, its decay proceeds via two gluon hadronization, which is similar 
to the second step of the Xc2 decay. Therefore, a pure 2~^~^ glueball may have flavor- 
symmetric decays. Furthermore, the 2~^~^ glueball, if lying in the 2230 MeV region, 
can only have little mixing with nearby L=3 2~^~^ quarkonium states, because these 
qq states have vanishing wave functions at the origin due to high angular momentum 
barrier which will prevent the qq pair from being annihilated into gluons and then 
mixed with the 2"*"+ glueball. This might explain why ^(2230) has flavor-symmetric 
couplings to TT+TT^ and K~^K~, if it is nearly a pure 2^~^ glueball. In addition, the 
gluon hadronization leads to many decay modes for Xco and Xc2, therefore the 2++ 
glueball may also have many decay modes. In comparison, the observed branching 
ratios for (^(2230) — > tt+tt", K~^K~, KsKs, pp may not exceed 6 percent. This 
might be very different from the conventional qq mesons, which usually have some 
dominant two-body decay modes. 

Above discussions indicate that the decay pattern of glueballs could be rather 
complicated, and a deeper theoretical understanding is needed to reduce the uncer- 
tainties. As for the glueball mass spectrum, despite the remarkable progress made in 
lattice QCD calculations [1U|^|^6|, uncertainties in estimating glueball masses are 
still not small. For instance, for 0^"*" glueball UKQCD group gives M = 1550±50MeV 
H, while IBM group gets M = 1740 ± 71MeV and F = 108 ± 29MeV [|| . 



Another progress in the lattice calculation is the glueball matrix elements. For ex- 
ample, a calculation for < 0\Tr{g'^G ^i^G fj_i^)\G > predicts a branching ratio of 5 x 10~^ 
in J/ip radiative decays for O"*""*" glueball |3^. This may provide useful information 
on distinguishing between /o(1500) and 6'(1720), or other possible candidates for 0"*"+ 
glueball. If /o(1500) is the 0"*"+ glueball, it should have some important decay modes 
(e.g., 47r) to show up in J/ip radiative decays. 

In summary, while the situation in searching for glueballs via charmonium decays 
is very encouraging, especially with the r-charm factory in the future, more theoretical 
work should be done to make more certain predictions on the glueball mass spectrum, 
the widths, the transition matrix elements, and in particular the decay patterns. 
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IV. Prompt Charmonium Production at Tevatron and Fragmentation of 
Quarks and Gluons 

The study of charmonium in high energy hadron collisions may provide an important 
testing ground for both perturbative QCD and nonperturbative QCD. 

In earlier calculations |Q, in hadronic colhsions the leading order processes 

99 99, 9Xc{Xc ^ 7^), <19 ^ QXciXc 7^), (8) 

were assumed to give dominant contributions to the cross section. But they could 
not reproduce the observed data for charmonium with large transverse momentum. 
This implies that some new production mechanisms should be important. These are 
the quark fragmentation and gluon fragmentation. 



1. Quark Fragmentation 

In essence the quark fragmentation was first numerically evaluated in a calculation 
for the decay ijjcc by Barger, Cheung, and Keung^^(for other earlier 



discussions on fragmentation mechanisms see ref. |^). This decay proceeds via 
cc, followed by the splitting c ^ ipc or c —>■ ipc (see Fig.l), of which the rate is two 
orders of magnitude larger than that for Z^ — > '?/'5f(7|^, because the fragmentation 
contribution is enhanced by a factor of {Mz / mc)'^ due to the fact that in fragmentation 
the charmonium (cc bound state) is produced with a seperation of order l/irtc rather 
than 1/mz as in the previous short-distance processes, e.g., Z° — > ijjgg. 



Fig. 1 The quark fragmentation mechanism. il) is produced by the charm quark 
sphtting c ipc. 

These numerical calculations, which are based on the fragmentation mechanisms, 
can be approximately (in the limit nic/mz 0) re-expressed in a more clear and 
concise manner in terms of the quark fragmentation functions, which were studied 
analytically by Chang and Chen[Q[^, and by Braaten, Cheung, and Yuan^^. The 
quark fragmentation functions can be calculated in QCD using the Feynman diagram 
shown in Fig. 1. For instance, the fragmentation function Dc^^{z, /i), which describes 
the probability of a charm quark to split into the J/ijj with longitudinal momentum 
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fraction z and at scale /i, is given byP^ 

where /i = 3mc and R(0) is the radial wave function at the origin of J/ip. Large 
logarithms of fi/rric for /i = 0{mz) appearing in Di^^{z,fi) can be summed up by 
solving the evolution equation 

i^-g-Di^A^,f^) = Y. J —Pi^ji^/yyf^)Dj^Ay^f^)y (10) 

where Pi^j{x, fi) is the Altarelli-Parisi function for the splitting of the parton of type 
i into a parton of type j with longitudinal momentum fraction x. The total rate for 
inclusive ip production is approximately 

r(Z° ^^ + X)= 2f(Z° ^ cc) dzDc^Jz, 3me). (11) 

Jo 

Then the branching ratio for the decay of into ip relative to decay into cc is 

^ 0.0234„,,2».,« . 2 X 10-, (12) 



which agrees with the complete leading order calculation of Z'^ ipcc in Ref . [HO 



Using the fragmentation functions, the production rates of the Be meson are 
predicted in Ref. ||2|. E.g., the branching ratio of Be in Z^ decay is about R ^ 
7.2 X 10-5 (see also Ref.|§). 

The quark fragmentation functions to P-wave mesons have also been calculated 

El 



2. Gluon Fragmentation 

As the quark fragmentation, the gluon fragmentation may also be the dominant 
production mechanism for heavy quark-antiquark bound states (e.g., charmonium) 
with large tansverse momentum. 

In previous calculations, e.g. in the gluon fusion process, charmonium states with 
large Pt were assumed to be produced by short distance mechanisms, i.e., the c 
and c are created with transverse seperations of order 1/Pt, as shown in Fig. 2 (a) for 
99 ~^ Vc9- However, in the gluon fragmentation mechamism the rjc is produced by the 
gluon splitting g rjc9 (while J /ip is produced by (? — > ''P99), as shown in Fig.2(b). 

For gluon fragmentation, in the kinematic region where the virtual gluon and ?7c 
are colinear, the propagator of this gluon is off shell only by an amount of order mc, 
and enhances the cross section by a factor of P^/m1. If Pt is large enough, this will 
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overcome the extra power of the couphng constant as, as compared with the short 
distance leading order process gg — ^ r]cg. The gluon fragmentation functions were 
calculated by Braaten and YuanPTf 



1 -\m\' 

7271 m' 



rf^D„(z,2m,) = _a,(2m,)2^_^, (13) 



' dzDg^^z, 2m,) = (1.2 X I0~^)asi2m,f^-^^^. (14) 



c 



where the latter is estimated to be smaller than the former by almost an order of 
magnitude. 



Fig. 2 (a) A Feynman diagram for gg ccg that contributes to rjc production 
at order a^; Fig. 2(b) A Feynman diagram for gg —^ ccgg that contributes 
to rjc production at order a^. For the virtual gluon at large Pr, with qq =- 
0{Pt), ~ O(m^), the contribution is dominant. 



The gluon fragmentation into P-wave heavy quarkonium was also studied ||48 
The P-wave state (e.g., Xcj) can arise from two sources i.e. the production of a color- 
singlet P-wave state, and the production of a cc pair in a color-octet S-wave state, 
which is then projected onto the Xcj wave functions. With two parameters which 
characterize the long-distance effects, i.e., the derivative of P-wave wavefunction at 
the origin and the probability for an S-wave color-octet cc pair in the color-singlet Xcj 
bound state, the fragmentation probabilities for a high transverse momentum gluon 
to split into Xco, Xa, Xc2 are estimated to be 0.4 x 10"'*, 1.8 x 10""^, 2.4 x 10~^, 
respectively [^8[|. They could be the main source of Xcj production at large Pt in pp 
colliders. 

Since fragmentating gluons are approximately transverse, their products are sig- 
nificantly polarized. Cho, Wise, and Trivedi |]50[] find that in gluon fragmentation to 



Xcj(l-P) followed by Xcj — ^ iJ/"^ the helicity levels of Xci, Xc2, and J/ip are populated 
according to certain ratios, e.g., -D^^° : D]^}^^ ^1:1, DjJ^ '■ D^jj^^ ~ 1 : 3.4. 

The gluon fragmentation to J^'" = 2^+ quarkonia was also studied, and 

these D-wave state's polarized fragmentation functions were computed . 
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3. The ifj' surplus problem at the Tevatron 

In 1994, theoretical calculations were compared with data on inclusive J /ip and ip' pro- 
duction at large transverse momentum at the Tevatron |5l|, where large production 



cross sections were observed. The calculations include both the conventional leading 



order mechanisms and the charm and gluon fragmentation contributions [p^ . 

For if) production both fragmentation directly into if) and fragmentation into Xc 
followed by the radiative decay Xc — ^ + 7 ^'^^ considered. Fragmentation functions 
for 5f — > '?/', c — '?/', g ^ Xc, c ^ Xc, and j ^ ip are used. 

These calculations indicate that 

(1) Fragmentation dominates over the leading-order mechanisms for Pt > 5 GeV. 

(2) The dominant production mechanism by an order of magnitude is gluon frag- 
mentation into Xc followed by Xc — *■ 7V'- 

For ip' production the fragmentaions g ^ ip', c ^ 7 — > ip' and the leading 
order mechanisms are included but no contribution from any higher charmonium 
states is taken into consideration. The dominant production mechanisms are gluon- 
gluon fusion for Pt < 5GeV, and charm quark fragmentaion into ip' for large Pt- 
However, the calculated production cross section of ip' is too small by more than an 
order of magnitude (roughly a factor of 30) This serious disagreement, the 

so-called ip' surplus problem, has caused many theoretical speculations. 

The radically excited 2^Pi 2 (xci(2-P) and Xc2{'^P)) states have been suggested to 
explain the xp' surplus problem |^ . These states can be produced via gluon 

and charm fragmentaion as well as the conventional gluon fusion mechanism, and 
then decay into 'yip' through El transitions. Large branching ratios of B{xcj{'^P) ^ 
ip'{2S) + 7) = (5 ~ 10)% (J=l, 2) are required to explain the ip' production 
enhancement. Within the potential model with linear confinment, the masses of 
these 2P states are predicted to be, e.g., M(xco(2P)) = 3920MeV, M(xci(2P)) = 
3950MeV, and M(xc2(2P)) = 3980Mel^ 0, therefore OZl-allowed hadronic decays 
like Xco(2P) ^ DD, Xci{2P) D*D + c.c, and Xc2(2P) ^ DD, D*D + ex., can 
occur. It is not clear whether these hadronic widths are narrow, making the brancing 
ratios -B(xcj(2P) — >■ ip'{2S) + 7) large enough to explain the ip' production data. 

One possibility is that since decays Xc2(2P) DD, D*D + c.c. proceed via L = 2 
partial waves, they could be suppressed |Q. These OZl-allowed hadronic decays are 
estimated in a fiux-tube model and they could be further suppressed (aside from 
the D-wave phase space for Xc2{2P)) due to the node structure in the radial wave 
functions of excited states With suitble parameters used, the widths of Xc2{2P) 



and Xci(2P) could be as narrow as F (1 ~ 10)MeV. 

There is another possibility that the Xci{2P) could lie bellow the D*D threshold, 
and then with roughly estimated r{xci{2P) — > light hadrons)~ 640 KeV, r{xci{2P) 
'yip') ^ 85 KeV, one could get B{xci{2P) -fip') ^ 12% i|. This possibility relies 
on the expectation that color screening effect of light quark pair on the heavy Q-Q 
potential, observed in lattice QCD calculations, would lead to a screened confinement 
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potential which makes the level spacings of excited cc states lower than that obtained 
using the linear potential (e.g., ip{4:160) and ip{A415) could be 4S and 5S rather than 
2D and 4S states, respectively). 

Moreover, it is suggested that the ccg hybrid states could make a significant con- 
tribution to J/ip and ip' signals at the Tevatron||S^, since the color octet production 



mechanism is expected to be important, and hybrid states contain gluonic excitations 
in which the cc are in the S-wave color octet configuration. In particular, the negative 
parity hybrid states, including (0, 1, 2)""'", lying in the range 4.2±0.2GeV, could 
be a copious source of J/ip and through radiative and hadronic transitions. 

4. Color Octet Fragmentation Mechanism 

Baseed on a general factorization analysis of the annihilation and production of heavy 



quarkonium ||6^[^, Braaten and Fleming proposed a new mechanism, i.e. the color- 



octet fragmentation for the J/ifj and ijj' production at large Pt \Q2 . 

In the framework of NRQCD theory ||61]|[0, which is based on a double power 
series expansion in the strong coupling constant as and the small velocity parameter 
V of heavy quark, the fragmentaion functions can be factored into short-distance 
coefficients and long-distance factors that contain all the nonperturbative dynamics 
of the formation of a bound state containing the cc pair. E.g., for g^ip fragmentation 

Dg^^z, fi)=Y^ dn{z, ^l)<Q\Ot\Q>, (15) 

n 

where O'^ are local four fermion (quark) operators. 

For the physical ip state, the wavefunction can be expressed as Fork state decom- 
positions which include dynamical gluons and color-octet {QQ)s components 

\ij> = OmQQU'S,)>+0{v)\{QQ)sCPj)9> 

+ Oiv')\iQQ)sCS^)gg >+■■■. (16) 

Therefore there are two mechanisms for gluon fragmentation into ip: 

(1) Color-singlet fragmentation g* -^ccgg. Here cc is produced in a color-singlet 
^Si state. The matrix element < Ofl^Si) > is of order mj?f^, which is related to 
the Fork state \{cc)i(^Si) > in ip, so the contribution to fragmentation function is of 
order a^v^. 

(2) Color-octet fragmentation g* cc. Here cc is produced in a color-octet ^5*1 
state. The matrix element < 0^{^Si) > is of order m^u'', which is related to the 
Fork state \{cc)s{^Si)gg > in ip, so the contribution to fragmentation function is of 
order a^f ^. 

It is clear that the color-octet fragmentation g* — > cc is enhanced by a factor of 
~ a~'^ from the short-distance coefficients, and suppressed by a factor of ~ f ^ from 
the long-distance matrix elements, as compared with the color-singlet fragmentation. 
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Since for charmonium ^ 0.25 ~ 0.30 is not very small, the color-octet fragmenta- 
tion could be dominant in some cases, e.g., in the ip' production at large transverse 
momentum. 

In the case of , if the observed large cross section is really due to color-octet 
fragmentation the matrix element < Og (^^i) > will be determined by fitting the 
CDF data on the production rate at large Pt 

< Ot\^Si) >= 0.0042G'eV^ (17) 

while the color-singlet matrix element < Of {^Si) > is determined by the ip' leptonic 
decay width which is related to the wave function at the origin 

< OfCSi) >^ — |i?^f = O.UGeV^ (18) 

The color-octet matrix element is smaller by a factor of 25 than the color-singlet 
matrix element, consistent with suppression by f^. Therefore the color-octet frag- 
mentation could be a possible solution to the ip' surplus problem. 

The color-octet fragmentation will also make a substantial contribution to the 
J/ijj production at large Pt, and may compete with gluon fragmentation into Xcj 
followed by Xcj iJ/i^- 

The color-octet fragmentation mechanism might be supported by the new data 
from CDF[^ and DO|Q at the Tevatron. New results for the fraction of J /ip which 



come from the radiative decay of Xc are (see Fig.3 for the CDF result) 

CDF : f'l'^ = (32.3 ± 2.0 ± 8.5)% (Py^'^ > AGeV, \r]^/^\ < 0.6), (19) 

DO : ///^ = (30 ± 7 ± 10)% [P^^^ > 8GeV, \r]-^/^\ < 0.6). (20) 



Fig.3 The fraction of Jf-ip from Xc as a function of P^^^' with the contribution 
from b quark's removed, measured by CDF |63]. 
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This implies that the majority of prompt J /ip at large Pt do not come from Xc-, and 
the gluon fragmentation into Xc is not the dominant mechanism for J/ip production 
at large Pt- The observed production cross section of J/tfj from Xc is in reasonable 
agreement with the theoretical calculations while the direct J/ip production cross 
section is a large factor above the prediction (see Fig.4 for the CDF result). 



Fig. 4 Differential cross sections of prompt J/'0 as a function of P^^^ with the 
contribution from b quark's removed, measured by CDF The dotted curve 
repesents the total fragmentation contribution (but without the color-octet 
fragmentation) , and the dashed curve represents the leading-order contribution 
||. 

Although this result might favor the color-octet fragmentation mechanism for the 
direct production of J/ip., it is still premature to claim that it is the real source of 
J/ip production. 

In order to further test the color-octet fragmentation mechanism for the produc- 
tion of J/ip and, in particular, ifj' at the Tevatron, some studies are required. First, 
the produced ifj' should be transversely polarized and the experimental obser- 



vation of a large transverse if)' spin alignment would provide strong support for the 
color-octet production mechanism of ip' . Another important test is to apply the same 
mechanism to the hb systems. The integrated and differential production cross sec- 
tions for the T(IS'), T(25'), T(35') have been measured by both CDF and DO 
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The production rates are generally found to be higher than that with color-singlet 
fragmentations. The color-octet production mechanism does help to explain some of 
the discrepancies [ |6^ . 

In this connection it is worthwhile to note that the problem of J /ip and especially 
ip' surplus production has also been observed by the fixed target experiment (e.g., in 
800 GeV proton-gold collisions) . In collisions at lower energies, fragmentation is 
not expected to be dominant. It is not clear whether the ip' surplus observed both at 
the Tevatron and fixed-target (with the same enhancement factor of about 25 relative 
to the expected production rates) has the same origin or not. Further experimental 
and theoretical investigations are needed. 



V. Some Results in Charmonium and Open Charm Physics 



Here some theoretical results on charmonium and open-charm hadrons are reported. 
• The QQ spin dependent potential. 

There have been many discussions about the QQ spin dependent potential (see e.g. 
miUlO)- A formula for the heavy- quark- ant iquark spin-dependent potential 



is given using the techniques developed in heavy-quark effective theory |]7T|. The lead- 
ing logarithmic quark mass terms emerging from the loop contributions are explicitly 
extracted and summed up. There is no renormalization scale ambiguity in this new 
formula. The spin-dependent potential in the new formula is expressed in terms of 
three independent color-electric and color- magnetic field correlation functions, and it 
includes both the Eichten-Feinberg formula and one-loop QCD result as 



special cases. For hyperfine splittings with ^jjg = 200 — 500MeV^, the new formula 
gives J7| M{J/iIj) - M{r]^) ^ 110 - UOMeV, M(T) - M(r/fe) ^ 45 - 50MeV, and 

MCPi) - M^Pj) ^ 2 - AMeV (21) 



for cc, which is larger than the present E760 result (~ O.dMeV) [|7^], and other 
theoretical predictions (e.g. [^). But this tiny mass difference may be sensitive to 
other effects, e.g., the coupled-channel mass shifts. 

A set of general relations between the spin-independent and spin- dependent poten- 
tials of heavy quark and antiquark interactions are derived from reparameterization 
invariance in the Heavy Quark Effective Theory [^. They are useful in understand- 
ing the spin-independent and spin-dependent relativistic corrections to the leading 
order nonrelativistic potential. 

• Relativistic corrections to QQ decay widths and the determination of as{mQ). 
Charmonium mass spectrum and decay rates can be very useful in determining the 
QCD coupling constant a^. In recent years remarkable progresses have been made 
in lattice calculations ||76[|7^. On the other hand, many decay processes may be 
subject to substantial relativistic corrections, making the determination of as quite 
uncertain 
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The decay rates of V 3g and V e^e~ for V = J /ip and T may be expressed 
in terms of the Bethe-Salpeter amphtudes, and to the first order relativistic correction 
and QCD radiative correction it is found that |8( 



V{y 3g) = — ^ d q[l- 2.95— ]^sch{q) U X 22 



where Sc = 3.7, Sh = 4.9 (defined in the MS scheme at the heavy quark mass 
scale) ^9] . This resuh shows exphcitly that the relativistic correction suppresses 
the gluonic decay much more severely than the leptonic decay. Using the meson 
wavefunctions obtained by solving the BS equation with a QCD-inspired interquark 



potential, and the experimental values of decay rates |81|, it is found that 



a,(mc) = 0.26 - 0.29, ^.(mf,) = 0.19 - 0.21, (23) 

at TTLc = 1.5 GeV and = 4.9 GeV. These values for the QCD coupling constant 
are substantially enhanced, as compared with the ones obtained without relativis- 
tic corrections. However, it should be emphasized that these numerical results can 
only serve as an improved estimate rather than a precise determination, due to large 
theoretical uncertainties related to the scheme dependence of QCD radiative correc- 



tions and higher order relativistic corrections. This result is consistent with that 



obtained using finite size vertex corrections [p^j . 
• Heavy meson decay constants. 

Discussions on the heavy meson decay constants are very extensive. In the frame- 
work of heavy quark effective theory (HQET), QCD sum rules are used to estimate 



the nonperturbative effects ||8J] [§6| |^ . The first systematic investigation was 
given in [^, and a further improvement was obtained by seperating the subleading 
order from the leading one Bl 



In a recent work |Q the SU(3) breaking effects in the leading and subleading 
parameters appeared in the heavy quark expansion of decay constans of the heavy- 
light mesons are systematically analyzed to two loops accuracy using QCD sum rules. 
It is found that the SU(3) breaking effects in the decay constant of the pseudoscalar 
are respectively 

fBjfB = 1.17 ± 0.03, foJfD = 1.13 ± 0.03. (24) 
These results are in agreement with recent lattice QCD calculations WTj . In addition. 



the ratios of vector to pseudoscalar meson decay constants are found to be 

= fBjfB = 1.05 ± 0.02, (25) 
and the SU(3) breaking effect in the mass is about 82 ± SMeV. 
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Another approach to estimating nonperturbative effects on heavy mesons is to 
combine HQET with chiral perturbation theory . In the framework of the heavy- 
hght chiral perturbation theory (HLCPT) the heavy meson decay constants are dis- 
cussed and the effects of excited states on the chiral loop corrections are further 
considered PT[] . 

In a recent work the vector meson contributions are introduced in HLCPT and 
the lagrangian and current to the order l/mg are constructed |92|. With this, to 
the order l/A^^^, {-^csb is the chiral symmetry breaking scale), corrections to fo and 
/b arising from coupled-channel effects to order l/rric and l/rrib are calculated. At 
the tree level in HLCPT, using the relativistic B-S equation with kernel containing 
a confinement term and a gluon exchange term in a covariant generalization of the 

when rriQ 



Coulomb gauge [Q, the decay constants /), 



(0) 



and fP 



oo as well as 

the l/rriQ corrections are calculated. HLCPT and the heavy quark effective theory 
(HQET) are matched at the scale Acsb- Adding the perturbative and nonperturbative 
contributions the values for fc and fs are found to be 



fo^fB^ 200 MeV, 



(26) 



which is in agreement with lattice calculations . 

We now turn to some new experimental results in open charm physics. The CLEO 
Collaboration has given following results. 

• More accurate or the first measurements of decays 



Channel 


B(%) 


PDG(%) 


K+K- 


0.455±0.029±0.032 


0.454±0.029 




0.048±0.012±0.013 


0.11±0.04 




0.074±0.010±0.018 


0.089±0.025 




<0.063 at 90 CL% 






0.107±0.030 





The theoretical prediction for B{K^K ) is in the range 0.14-0.6, and for 
is 0-0.3. 

• Observation of the Cabibbo suppressed charmed baryon decay of A^ 
pK^K~ , compared with Ajf pK~7c^ [^. 

B{p<p)/B{pKTx) = 0.024 ± 0.006 ± 0.003, 

B{pKK)/B{pK7r) = 0.039 ± 0.009 ± 0.007, 
B{p(f)) /B{pKK) = 0.62 ± 0.20 ± 0.12. 
The theoretical predictions range from 0.01 to 0.05 for B{p(f))/B{pK'n') 



BiK'^K^) 
pep and 

(27) 

(28) 
(29) 



Measurement of the isospin-violating decay D 



0.062l°:°20 ^ Q_Q22. 



(30) 
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This isospin-violating decay is expected to proceed through OZI-allowed decay 
D*^ Dff] (via the ss component in 77) and the — vr^ mixing This decay also 
imphes that D*~^ has natural spin-parity (most likely 1^). 

• Measurement of the relative branching ratios of to rje^u and rj'e^u, com- 
pared to 



B{D+ rje+u) 
B{D+ (f)e+u) 

B{D+ T]'e+iy) 



1.24 ±0.12 ±0.15. 



0.43 ±0.11 ±0.07. 



These results favor the prediction of the ISGW2 model |99 



Measurement of decay branching ratios relative to 



(31) 
(32) 
E00|. 



Decay Mode 


events 


B/BiE+ E'n+n+) 




119±23 


1.18±0.26±0.17 




61±17 


0.92±0.27±0.14 




61±15 


0.58±0.16±0.07 


e-7r+7r+ 


131±14 


1.0 



There are also some experimental results from the ARGUS Collaboration. 
• Leptonic branching ratios of [101]. 

B{D° e+z/eX) = 6.9 ± 0.3 ± 0.5%, 

B{D^ li^u^X) = 6.0 ± 0.7 ± 1.2%. 



(33) 
(34) 



These values are smaller than the world average values |81[ . 

• Measurement of the decay D^(2573) —>■ D^K^ ||102|| . The observed mass and 
width r = (10.4 ±8.3 ±3.0) MeV of this resonance are consistent with that obtained 
by CLEO. 

• Evidence for the A*+(2593) production p3] . 

Finally, BES Collaboration has reported the leptonic branching ratio of Ds using 

(35) 



;i48 ± 18 ± 13) D, events [|T04| 

B{Dt ^ e+i/eX) 



4.6-1.2J/0- 



VI. Conclusions 



While impressive progress in experiment has been made in physics in the charm 
energy region, some theoretical issues need to be clarified. 

The new data and puzzles in exclusive hadronic decays of J/tp and ip' give new 
challenges to the theory of hadronic decays. 
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With the new observation for ^(2230) and /o(1500), the situation in searching 
for glueballs is encouraging, but theoretical uncertainties related to the properties of 
glueballs still remain and need to be further reduced. 

For the prompt production of charmonium at large transverse momentum, gluon 
and quark fragmentations dominate over leading-order parton fusions. Color-singlet 
fragmentation is not the dominant mechanism for J/%1) and if)' production. Color- 
octet fragmentation seems to be important to explain the J/%1) and, in particular, the 
ip' excess, but further tests are required. The mechanism of charmonium production 
at fixed-target also needs studying. 

The study of open charm physics is in continuous progress. This is important for 
testing the Standard Model and understanding both perturbative and nonpcrturba- 
tive QCD. 

In the future, with new experiments at e'^e~ colliders, hadronic colliders, fixed 
target, and, in particular, at the proposed r-charm factory, and with the theoretical 
progress in lattice QCD and other nonperturbative methods, a deeper understanding 
of physics in the charm energy region will be achieved. 

Acknowledgements 

I would like to thank my colleagues, in particular, Y.B. Dai, Z.X. He, T. Huang, Y.P. 
Kuang, J.M. Wu, and H. Yu for many very helpful discussions and suggestions. I 
would also like to thank Y.F. Gu, S. Jin, J. Li, W.G. Li, C.C. Zhang, Z.P. Zheng, 
and Y.C. Zhu for useful discussions on experimental results. Thanks are also due to 
H.W. Huang, J.F. Liu, Y. Luo, and especially C.F. Qiao for their help in preparing 
this report. I also wish to thank V. Paradimitriou for providing me with the new 
data from CDF. 

References 

[1] Mark II Collaboration, M.E.B. Franklin et al, Phys. Rev. Lett. 51 (1983) 963; 
Mark III CoUabotation, W.H. Toki, SLAC-PUB-5093 (1989). 

[2] J. Li, EES Collaboration, these proceedings. 

[3] EES Collaboration, J.Z. Eai et al, High Energy Physics and Nuclear Physics, 
19 (1995) 93. 

[4] S.J. Erodsky and CP. Lepage, Phys. Rev. D24 (1981) 2848. 
[5] P. Freund and Y. Nambu, Phys. Rev. Lett. 34 (1975) 1645. 
[6] W.S. Hou and A. Soni, Phys. Rev. Lett. 50 (1983) 569. 



19 



[7] S.J. Brodsky, G.P. Lepage and S.F. Tuan, Phys. Rev. Lett. 59 (1987) 621. 
[8] S.S. Pinsky, Phys. Lett. B236 (1990) 479. 

[9] M. Chaichian and N.A. Tornqvist, Nucl. Phys. B323 (1989) 75. 

[10] G. Karl and W. Roberts, Phys. Lett. B144 (1984) 263; G. Karl and S.F. Tuan, 
Phys. Rev. D34 (1986) 1629. 

[11] Crystal Ball Collaboration, R. Lee, Ph.D. Thesis, Stanford University, SLAG 
Report-282 (1985). 

[12] Y.F. Gu and S.F. Tuan, Mod. Phys. Lett. AlO (1995) 615. 

[13] M. Anselmino et al, Phys. Rev. D44 (1991) 1597; D50 (1994) 595. 

[14] K.T. Ghao, Y.F. Gu and S.F. Tuan, to appear in Gommun. Theor. Phys. 

[15] S. Jin, BES Collaboration, Proceedings of the Workshop on Hadron Physics at 
e+e" CoUider, Beijing, 1994, edited by: X.L. Fan, W.G. Li and H.L. Ni, p.77. 

[16] R.M. Baltrusaitis et al, Phys. Rev. Lett. 56 (1986) 107. 

[17] P.D. Barnes et al, Phys. Lett. B309 (1993) 469. 

[18] K.T. Ghao, Gommun. Theor. Phys. 24 (1995) 373 (|hep-ph/ 950240^ ). 

[19] T. Huang, S. Jin, D.H. Zhang and K.T. Ghao, BIHEP-TH-95-11, PUTP-95-09. 

[20] S. Goldfrey, R. Kokoski and N. Isgur, Phys. Lett. B141 (1984) 439. 

[21] H. G. Blundell and S. Godfrey, OGIP/G 95-11, |hi^i/9508264 

[22] G.S. Bah et al, Phys. Lett. B309 (1993) 378. 

[23] Z. Bai et al, Phys. Rev. Lett. 65 (1990) 2057. 

[24] J.E. Augustin et al, Phys. Rev. D46 (1992) 1951. 

[25] A.M. Ma, BES Collaboration, Ph.D. Thesis, 1995. 

[26] M. Chanowitz, Proceedings of GCAST (World Laboratory) Symposium on 
Charm Physics, Beijing, 1987, edited by: M.H. Ye and T. Huang (Gordon and 
Breach, 1988) p.l61. 

[27] G.J. Gounaris and H. Neufeld, Phys. Lett. B213 (1988) 541; G.J. Gounaris and 
J.E. Paschalis, Phys. Lett. B251 (1990) 634. 



20 



[28] C. Amsler et al, Phys. Lett. B340 (1994) 259; B342 (1995) 433; B291 (1992) 
347. 

[29] S. Abatzis et al, Phys. Lett. B324 (1994) 509. 
[30] D. Aide et al, Phys. Lett. B201 (1988) 160. 

[31] D.V. Bugg, Proceedings of the Second International Symposium on Medium 
Energy Physics, Beijing, China, 1994, edited by W.Q. Chao and P.N. Shen 
(World Scientific, Singapore, 1995)p.51. 

[32] S.S. Gershtein et al, Z. Phys. C24 (1984) 305. 

[33] C. Amsler and F.E. Close, Phys. Lett. B353 (1995) 385. 

[34] N. T5rnqvist, Z. Phys. C61 (1994) 525. 

[35] C. Michael and M. Teper, Nucl. Phys. B314 (1989) 347. 

[36] H. Chen, J. Sexton, A. Vaccarino and D. Weingarten, Nucl. Phys. B(Suppl.) 34 
(1994) 357; D. Weingarten, talk given at the International Conference on Lattice' 
95. 



[37] Y. Liang, K.F. Liu, B.A. Li, S.J. Dong and K. Ishikawa, UK/92-05, f^- 
I lat/ 930401^ 



[38] R. Baier and R. Riickl, Z. Phys. C19 (1983) 251; F. Halzen et al, Phys. Rev. 
D30 (1984) 700. 

[39] B. Guberina, J.H. Kwhn, R.D. Peccei, and R. Riickl, Nucl. Phys. B174 (1980) 
317. 

[40] V. Barger, K. Cheung, and W.Y. Keung, Phys. Rev. D41 (1990) 1541. 

[41] For the gluon fragmentation mechanism of ip production in decay, see K. 
Hagiwara, A.D. Martin, and W.J. Stirling, Phys. Lett. B267 (1991) 527; 
For earlier discussions on quark fragmentation see also: C.R. Ji and F. Amiri, 
Phys. Rev. D35 (1987) 3318; J.H. Kwhn and H. Schneider, Phys. Rev. D24 (1981) 
2996; Z. Phys. Cll (1981) 263. 

[42] C.H. Chang and Y.Q. Chen, Phys. Lett. B284 (1992) 127; Phys. Rev. D46 (1992) 
3845. 

[43] Y.Q. Chen, Phys. Rev. D48 (1993) 5181. 

[44] E. Braaten, K. Cheung, and T.C. Yuan, Phys. Rev. D48 (1993) 4230; Phys. Rev. 
D48 (1993) R5049. 



21 



[45] C.H. Chang, Y.Q. Chen, CP. Han, and H.T. Jiang, AS-ITP-94-24, f^- 
I ph/950740^ . 



[46 
[47; 
[48 
[49 
[50 
[51 

[52 
[53 

[54; 
[55 
[56 

[57; 

[58 
[59 
[60 
[61 
[62 
[63 

[64 
[65 
[66 
[67 



T.C. Yuan, Phys. Rev. D50 (1994) 5664. 

E. Braaten and T.C. Yuan, Phys. Rev. Lett. 71 (1993) 1673. 

E. Braaten and T.C. Yuan, Phys. Rev. D50 (1994) 3176. 

P.Cho and M.B. Wise, Phys. Lett. B346 (1995) 129. 

P.Cho, M.B. Wise, and S.P. Trivedi, Phys. Rev. D51 (1995) R2039. 

CDF Collaboration, F. Abe et al, Phys. Rev. Lett. 69 (1992) 3704; CDF Col- 
laboration, T. Daniels, Fermilab-Conf-94/136-E. 

M. Cacciari and M. Greco, Phys. Rev. Lett. 73 (1994) 1586. 

E. Braaten, M.A. Doncheski, S. Fleming, and M.L. Mangano, Phys. Lett. B333 
(1994) 548. 

D.P. Roy and K. Sridhar, Phys. Lett. B339 (1994) 141. 

F. E. Close, Phys. Lett. B342 (1995) 369. 

D. P. Roy and K. Sridhar, Phys. Lett. B345 (1995) 537. 
S. Godfrey and N. Isgur, Phys. Rev. D32 (1985) 189. 
P.R. Page, OUTP-94-34P, |hep-ph/9502204| . 

K.T. Chao and C.F. Qiao, PUTP-95-17. 

G. T. Bodwin, E. Braaten, and G.P. Lepage, Phys. Rev. D51 (1995) 1125. 
W.E. Caswell and G.P. Lepage, Phys. Lett. B167 (1986) 437. 

E. Braaten and S. Fleming, Phys. Rev. Lett. 74 (1995) 3327. 
CDF Collaboration, Fermilab-Conf-95/226-E. 

DO Collabotation, L. Markosky, Fermilab-Conf-95/137-E. 
CDF Collaboration, Fermilab-Conf-95/227-E. 
P. Cho and A.K. Leibovich, CALT-68-1988, |hep-ph/ 9505329 . 
M.H. Schub et al, Phys. Rev. D52 (1995) 1307. 



22 



[68] E. Eichten and F. Feinberg, Phys. Rev. D23 (1982) 2724. 

[69] D. Gromes, Z. Phys. C26 (1984) 401. 

[70] W. Buchmiillcr, Y.J. Ng, and S.H.H. Tye, Phys. Rev. D24 (1981) 3003; J. Pan- 
taleone, S.H.H. Tye, and Y.J. Ng, Phys. Rev. D33 (1986) 777; S.N. Gupta, S.F. 
Radford and W.W. Repko, Phys. Rev. D26 (1982) 3305. 

[71] Y.Q. Chen, Y.P. Kuang, and R.J. Oakes, Phys. Rev. D52 (1995) 264. 

[72] Y.Q. Chen and R.J. Oakes, NUHEPTH-95-05. 

[73] T.A. Armstrong et al, Phys. Rev. Lett. 69 (1992) 2337. 

[74] F. Halzen et al, Phys. Rev. D47 (1993) 3013. 

[75] Y.Q. Chen and Y.P. Kuang, Z. Phys. C67 (1995) 627. 

[76] A. EL Khadra, G. Hockney, A. Kronfeld, and P. Mackenzie, Phys. Rev. Lett. 69 
(1992) 729; G.P. Lepage and J. Sloan, NucL Phys. B (Proc. SuppL) 34 (1994) 
417; C.T.H. Davies et al, Phys. Rev. D50 (1994) 6963; C.T.H. Davies et al, 
OHCTPY-HEP-T-94-013. 

[77] For a recent review on lattice QCD theory, see: J. Shigemitsu, Proceedings of 
the XXVII International Conference on High Energy Physics, Glasgow, UK, July 
1994, p.135. 

[78] P. B. Mackenzie and G. P. Lepage, Phys. Rev. lett. 47 (1981) 1244. 

[79] W. Kwong, P. B. Mackenzie, R. Rosenfeld, and J. L. Rosner, Phys. Rev. D37 
(1988) 3210. 

[80] K.T. Chao et al, Proceedings of the International Conference on Quark Con- 
finement and the Hadron Spectrum, Como, Italy, June 1994, p. 306; K.T. Chao, 
H.W. Huang, and Y.Q. Liu, Phys. Rev. D, to appear. 

[81] Particle Data Group, L. Montanet et al, Phys. Rev. D50 (3-1) (1994) 1171. 

[82] S. Brodsky G. P. Lepage, and P. Mackenzie, Phys. Rev. D28 (1983) 228. 

[83] H.C. Chiang, J. Hiifner, and H.J. Pirner, Phys. Lett. B324 (1994) 482. 

[84] M. Neubert, Phys. Rev. D46 (1992) 1076. 

[85] V. Elesky and E. Shuryak, Phys. Lett. B276 (1992)191; E. Bagan et al, Phys. 
Lett. B278 (1992) 457. 

[86] P. Ball, Nucl. Phys. B421 (1994) 593. 



23 



[87] T. Huang and C.W. Luo, Phys. Rev. D50 (1994) 5775. 
[88] T. Huang and C.W. Luo, BIHEP-TH-94-39. 

[89] M.B. Wise, Phys. Rev. D45 (1992) 2188; T.M. Yan, H.Y. Cheng, C.Y. Cheung, 
G.L. Lin, Y.C. Lin, and H.L. Yu, Phys. Rev. D46 (1992) 1148; G.Burdman and 
J. Donoghue, Phys. Lett. B280 (1992) 287. 

[90] B. Grinstein et ai, NucL Phys. B380 (1992) 369. 

[91] A.F. Falk, Phys. Lett. B305 (1993) 268. 

[92] Y.B. Dai, X.H. Guo, C.S. Huang, and H.Y. Jin, AS-ITP-94-10. 

[93] H.Y. Jin, C.S. Huang, and Y.B. Dai, Z. Phys. C56 (1992) 707; Y.B. Dai, C.S. 
Huang, and H.Y. Jin, Z. Phys. C60 (1993) 527; C65 (1995) 87. 

[94] CLEO CoUaboration, CLEO CONE 95-20. 

[95] CLEO CoUaboration, CLEO CONE 95-16. 

[96] CLEO Collaboration, CLEO CONE 95-17. 

[97] P. Cho and M. Wise, Phys. Rev. D49 (1994) 6228. 

[98] CLEO CoUaboration, CLEO CONE 95-21. 

[99] D. Scora and N. Isgur, CEBAE-TH-94-14. 

[100] CLEO Collaboration, CLEO CONE 95-15. 

[101] ARGUS Collaboration, contributions for LP-95. 

[102] ARGUS Collaboration, contributions for LP-95. 

[103] ARGUS Collaboration, contributions for LP-95. 

[104] BES CoUaboration, contributions for LP-95. 



24 



